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ABSTRACT 

In this paper, we present a detailed high-resolution spectroscopic study of post main 
sequence stars in the Globular Gluster M68. Our sample, which covers a range of 4000 
K in Teff, and 3.5 dex in log(( 7 ), is comprised of members from the red giant, red hori¬ 
zontal, and blue horizontal branch, making this the first high-resolution globular cluster 
study covering such a large evolutionary and parameter space. Initially, atmospheric 
parameters were determined using photometric as well as spectroscopic methods, both 
of which resulted in unphysical and unexpected Tgff, log(( 7 ), ^t, and [Fe/H] combinations. 
We therefore developed a hybrid approach that addresses most of these problems, and 
yields atmospheric parameters that agree well with other measurements in the litera¬ 
ture. Furthermore, our derived stellar metallicities are consistent across all evolutionary 
stages, with ([Fe/H]) = —2.42 (cr = 0.14) from 25 stars. Chemical abundances obtained 
using our methodology also agree with previous studies and bear all the hallmarks of 
globular clusters, such as a Na-0 anti-correlation, constant Ca abundances, and mild 
r-process enrichment. 

Subject headings: stars: abundances — stars: evolution — stars: horizontal-branch - 
nucleosynthesis 


1. INTRODUCTION 


Along with M15 and M92 , M68 is one of the lowest metallicity Galactic globular clusters (GC), 
[Fe/H] ~ -2. j^ lHarrisI (llOodl. It has therefo re been included in many large-sample light-element 
abundance studies ie.g.. lBellazzini et al.ll2012l ). but has thus far been subjected to very few detailed 
chemical composition investigations. In an analysis of red horizontal branch (RHB) stars in the 
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metal-poor GC M15, iPreston et al.1 (|2006l ) found a declining difference between surface gravities 
determined from photometry and LTE spectrum analysis with increasing effective temperature in 
the range 5300 K < Tgff < 6300 K, a temperat ure range which embraces almost the entire RHB 
of that cluster. Contemporaneously, iLee et al.l (2005, hereafter Lee05) found larger differences 
in surface gravities among the cooler (Tgff ~ 4200 K) red giant branch (RGB) stars of M68. In 
combination, these results suggest that non-LTE (NLTE) over-ionization of neutral metals produces 
systematic errors in abundance analyses of cool, metal-poor red giants in globular clusters, and that 
abundances derived from RHB stars may provide a more accurate abundance scale for metal-poor 
stars of globular clusters and the Galactic halo. 

We decided to pursue this possibility by conducting a n expanded in vestigation of post main 
sequence stars in M68, similar to that for M15 reported bv ISobeck et al.l (2011, hereafter Sobll). 
Erom a purely observational point of view, M68 has many adv antages over M15. First, the RR 
Lyrae stars in M68 are ~0.2 mag. brighter than those in M15 (jWalkeii (1994, hereafter Wal94)). 
Second, all of our observations were done at the Las Campanas Observatory (LCO), and since M68 
has 5 = —27°, it transits the meridian near the zenith, whereas M15, with 5 = -|-12°, lies low in 
LCO’s northern sky. See Tabled] for more fundamental parameters of M68. 

In addition to being more accessible observationally, M68 differs from M15 in several other 
important respects. M68 has a relatively sparsely populated RHB comp ared to M15, which pos¬ 
sesses an extended blue horizontal branch (EHB) (jPurrell fc Harrislll993l ). This is also reflected in 
the horizontal branch ratio, defined to be the ratio of horizontal branch stars to red giant stars 


for b oth clusters. For M15, this ratio is 0.67, while for M68 it’s only 0.17 (|Zinn 


199611. Such differe nces in HB morphology have been studied extensively (e.g. iLee et al 


Harris 


1994 and 


Caputo et al.lll980h and attributed to variations in cluster age a nd metallicity, as well as stellar 
helium abundances and rotation (see Figs. 2-5 in lLee et al.lll994^ . Since M15 and M68 have very 
similar ages and metallicities, the differences in HB morphology could be attributed to disparate 
He contents. 


M68 and M15 also seem to exhibit systematically different abundance patterns. Lee05 found 
([Si/Fe])M68 ~ 0.6, while ([Si/Fe])Mi 5 '•^ 0-2 (Sobll). Si ab undances in ‘normal’ Pop H stars are 
equivalent to those of M68 RGBs fLee05. ICavrel et al.ll2004^ . Furthermore, M68 seems underabun- 
dant in Ti (Lee05), whereas overabundances of neutron capture elements, which vary from star 
to star, are f ound in the RHB and RGB stars of M15. Finally, M15 contains dusty red giants 
(|Bover et al.l (120061 )). surprising in view of the low metallicity, and interesting as unusually large 
mass loss during po st main sequence evolution has been advanced as an explanation for the EHB 
( D’Gruz et al. 1996 ). 


In this paper, we derive atmospheric parameters in a self-consistent fashion for RGB, RHB, and 
blue horizontal branch (BHB) stars in M68, which span about 4000 K in Tgif and 3.5 dex in log( 5 ). 
Figured] which shows a small spectral region for all evolutionary stages in our sample, highlights 
the difficulties of a self-consistent analysis over such a large parameter space. Several challenges, 
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such as the breakdown of photometric temperature calibrations, as well as the unphysicality of 
certain spectroscopic methodology assumptions, have to be overcome. After exploring photometric 
and spectroscopic methods to determine the values of Teg, log( 5 ), and [Fe/H], we develop a 
hybrid atmospheric analysis strategy that appears to yield reasonable parameters for the RGB, 
RHB, and BHB stars. Using these parameters, we present a detailed abundance analysis, which 
will allow us to gain insight about the chemical evolution of M68. 


OBSERVATIONS AND REDUCTIONS 


We obtained high resolution spectra of 11 red giant branch, 9 red horizontal branch, and 5 blue 
horizontal branch members of M68. All of our program stars were selected from the photometric 
survey of Wal94, whose V and B — V values are listed in in Table[2l Figure [2] shows the Walker color- 
magnitude diagram for M68, using symbol shapes and colors to distinguish the RGB, RHB, and 
BHB stars observed in this study. The observed RGB members were selected to represent stars from 
the giant branch tip down to the luminosity of the horizontal branch, while RHB candidates were 
chosen to cover most of their evolutionary stage up to the red edge of the RR Lyr gap (H—U ~ 0.45). 
BHB candidates were chosen to be between the blue RR Lyr edge (R — U ~ 0.25) and the 
domain [B — U ^ 0.2) in which stellar atmosnheric effects beein to distort the observed chemical 


Khalack et al. 

20] 

lA 

Behr 


2 OO 3 I and ref erences therein). Although M68 


has more than 40 known RR Lyr stars fe.g.. iGastellani et al.ll2003l and references therein), none of 
these were included in our study since spectrograph integration times would have been too long to 
acquire adequate data for these rapidly changing stars. 

Our spectra were gathered with the Magellan Inamori Kyocera Echelle (MIKE) spectrograph 


( Bernstein et al.ll2003l rl of the LCO Magellan Glay 6.5 m telescope. The spectrograph was config¬ 
ured with a 0.7" entrance aperture that yielded an ultimate resolving power of R = A/AA ~ 40,000 
for both the blue and red arms of instrument. The useful spectral coverage of the blue arm was 
3500—5000 A, and that of the red arm was 5000—9000 A. Scattered light and sky subtraction, as 
well as cosmic-ray filtering and flat field division were performed with software developed by S.A. 
Shectman (2004, unpublished). Wavelength calibrations were based on co-added hollow cathode 
Th-Ar spectra, obtained before and after each observation. Other one-dimensional extractions were 
completed using the apall package of IRAF@. 


^ http: //www.lco.cl/telescopes-information/magellan/instruments/mike 

®IRAF is distributed by the NOAO, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foundation. 
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3. LINE LISTS AND EQUIVALENT WIDTHS 
3.1. Atomic Line List 


Table [3] lists the lines of atomic species and their associated log{gf) values employed in this 
work. Measured equivalent width (EW) values, as well as the sou rces of our q/-values can be found 
in last column of this table. We call attention to Fe II, for which iMelendez &: Barbuvi ((20091) have 
proposed a renormalization of its transition probabilities based on some laboratory gf’s and an 
inverted solar spectrum analysis. These renormalized gf values are in general about 0.1 dex lower 
tha n those in the NIST A tomic Spectra Databas^, which were employed by us. If we had adopted 
the IMelendez &: Barbuvi results here, our Fe Il-based abundances would be larger by about 0.1 dex, 
which, in turn, would have decreased derived gravities by about 0.2 dex. We will return to this 
point later. 


3.2. Equivalent Width Measurements 


The equivalent width (EW) meas urements were done i n a semi-automated man ner using an 


IDL c ode (EW.pro) initially described in iRoederer et al.l ((20101) and further developed bv iBrugamver et al 


( 201 il l. This code allows the user to visually inspect either a Gaussian or Voigt minimization 
fit for each line, ensuring that any obviously blended or any otherwise undesirable line will not 
be measured. Additionally, the user can adjust the continuum which reduces any error possibly 
introduced by faulty normalization of the spectra. Several lines were picked ou t at random from 
all stars and re-measured using the SPECTRE code ( Fitzpatrick fc Sneden 1987 '^. Four members 
of our sample (RGB stars 256, 472 and RHB stars 403, 458) were analyzed entirely using SPEC¬ 
TRE. The EWs obtained in this fashion were then compared to those derived using the IDL code. 
On average, the difference between the EW values derived using EW.pro and those derived using 
SPECTRE are Aew = —1-18 {a = 3.18) mA. We therefore regard the differences in measured EW 
as negligible. 


The S/N, which is a function of wavelength and Teg of our programs stars, directly influenced 
the EW limitations for lines used in our analysis. A set of final S/N estimates for all of our stars 
can be found in Tabled! All S/N values given in this Table were estimated using the no routine of 
SPECTRE at around 6600 A in each star. For most stars, the lower EW cutoff was ~ 10 mA, while 
the upper limit was set at ~ 150 mA, depending on evolutionary state and line species. These 
limits, especially the lower cutoff, were not applied to star 334, as its S/N is about 70% lower than 
the average of our sample. For this star, any line that could be measured with a healthy degree of 


^ NIST is the U.S. National Institute of Standards and Technology; for the atomic line database see 
http://physics.nist.gov/PhysRefData/ASD/lines_form.html 

® Available at http://www.as.utexas.edu/~' chris/spectre.html 
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certainty, except the obviously saturated ones, was used. 


4. MODEL ATMOSPHERIC PARAMETERS 


4.1. Initial Parameters 


General information about M68, relevant to determining atmospheric parameters of our stars, 
can be found in Tabled! Special attention is called to the assumed distance modulus, (m — M)v, 
and the reddening, E{B — V), as they are of importance in determining photometric Tgfj and log( 5 ) 
values. 

Initial atmospheric parameters were determined from BVI photometry obtained by Wal94 
(see Table ED. To convert th e given fluxes to Tpff v alues, the (B-V) and (V-Ic) IRFM (infrared 
flux method) calibrations o f iRamfrez fc Melendea (120051 ! were used. Although there is a more 


recent calibration available (jCasagrande et al.l 120101 ) , only IRamfrez fc Melendea (|2005l ) include a 
calibration for giant stars. Figure [3| compares the Teg values obtained from (B-V) and (V-Ic) 
colors. In the RGB, these two indices give approximately the same answer, while they start to 
diverge in the RHB and BHB. This behavior can be explained by the lack of calibration stars at 
(V-Ic) < 0.6. Additionally, the difference between the V and Ic fluxes becomes insensitive to 
changes in temperature at (V-Ic) < 0.65, since the bandpasses of the respective filters are now in 
the temperature insensitive tail of the blackbody distribution. For these reasons, we chose (B-V) 
to be the sole color index in determining photometric temperatures. 

Initial surface gravity values were derived using the standard formula: 

log(5)* = 0.4(Mv* + BC - MboI©) + log(5)0 + 41og + log (^) ' 

For the solar values, MboI© = 4.75, Teg© = 5777 K, and log(g)© = 4.44 km s“^ were assumed. The 
assumed stellar mass in the above equation was estimated to be M* ~ 0.7 Mq, a result obtained 
from a isochrone calculation with M68 metallicity and age values (see 114.51 for more details). Since 
log( 5 )^ varies linearly with log(M*), an accurate value of stellar mass was not needed in this 
calculat ion. The bolometr ic correction (BC) for each star was calculated from the calibration 
given in I Alonso et al.l (|l999l i. This calibration, ho wever, does no t hold for stars with Tgfj > 6300 K. 
For stars exceeding this temperature. Figure 3 of Flower (jl996l ) was used to obtain an estimate of 


the stellar BC. Teg and log(g) values obtained in this fashion will be referred to as PHOT in texts 
and figures throughout the rest of this paper. 

Initial microturbulent velocities (^t) were esti mated to be 1.2 km s~ ^ for all RGB stars in our 
sample. For RHB stars, calibrations provided in iGratton et al.l (|l996l ) were used. Initial BHB 
esti mates were obtained by calculating a temperature dependent linear fit of a previous BHB 
study (jFor fc Snedenll2010l ). and applying the resulting calibration to our stars. Lastly, an initial 
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atmospheric metallicity of [Fe/H] 


—2.23 (from the 2010 editioij^ of [Harris 1996) was assumed. 


4.2. Final Parameters 


We employed three different methodologies to converge on our final set of parameters. In 
addition of the purely photometric PHOT parameter set defined above, we also derived purely 
spectroscopic parameters. These will be referred to as SPEC. Our final adopted set of parameters, 
which are a combination of photometry and spectroscopy, are called COMB. These designations 
will be used throughout the rest of this paper, including all plots. Abun dances for the various 
approaches were calculated using the latest version of MOOG ( Sneden 1973 except in the case of 
the RGB stars. In their temperature/gravity/metallicity regime, the major electron donors for the 
H~ continuous opacity are Fe and a elements. Since these elements are very dehcient in metal poor 
stars such as M68 RGB members, the H~ opacity decreases significantly, and scattering processes 
become important in the blue-UV spectral regions. Therefore, we employed a MOOG version 
incorporating Thompson scattering (see Sobll for more details) for this subgroup of our sample 


[X/H] abundances were calculated using the solar abundance recommendations of lAsphmd et al 
(|20n9l ). Our model at mospheres were interpol ated from ATLAS9 a-enhanced opacity distribution 
function model grids (ICastelli fc Knrucal2nr)3l ) using software developed by Andy McWilliams and 
Inese Ivans. 


To obtain the final atmospheric parameters for our stars, we decided to employ the analytical 
tools of the ‘classical’ spectroscopic approach, but deviated somewhat from its exact prescriptions. 
Specifically, we adopted photometric temperatures and then used the well known plots of individual 
Fe I line abundances as a function of excitation potential (EP) and as a function of reduced 
equivalent width, log(RW) = log(EW/A) to estimate microturbulent velocities. These photometric 
Teflf values were generally higher than the spectroscopic ones (see ^4.31 for more details), and an 
undesirable positive trend in the EP plot described above was obtained. This trend was mitigated 
by adjusting until acceptable correlations for both the EP and log(RW) plots were reached. 
Positive/negative changes in generally decrease/increase the abundances derived from strong 
lines; weak lines are generally not affected by changes in 

The log(( 7 ) values of our final approach were obtained by requiring equality between the abun¬ 
dances of [Fe l/H] and [Fe II/H]. This parameter determination process invariably changed the 
stellar metallicity, and therefore also implied photometric temperatures of the individual stars, 
forcing us to repeat the above process with these altered photometric values. Metallicities 
resulting from this second iteration differed very little from those of the first iteration, eliminating 
the need for a third iteration. Lastly, the metallicity of our atmosphere was equalled to ([Fe/H]) 


® http://www.physics.mcmaster.ca/' harris/mwgc.dat 
^ Available at http://www.as.utexas.edu/~ chris/moog.html 
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of each star. Atmospheric parameters obtained from the above methodology can be found in Ta¬ 
ble 0] under the COMB heading. Unless indicated otherwise, any further mention of atmospheric 
parameters will refer to these COMB values. 

Before highlighting two main successes of this approach, we review the methodology of the 
‘classical’ spectroscopic approach since we compare results of our final approach to those of the 
spectroscopic method below. Spectroscopic parameters (designated SPEC) were derived by: (1) 
requiring no trend with solar normalized abundances of Fe I and Fe II with excitation potentials 
(EPs) of different lines (giving Teg); (2) forcing ionization balance between Fe I and Fe II (giving 
log( 5 )); (3) demanding a correlation coefficient smaller than 0.01 between solar normalized line 
abundances and the logarithm of the reduced equivalent widths (log(RW)) (giving and (4) 
setting the atmospheric [Fe/H] equal to the resulting average normalized stellar iron abundance. 

The first motivation for our COMB approach is revealed in Figure 01 where we show trends 
of Fe and Ti abundances with Tgg. By design, the abundances of the neutral and ionized lines Fe 
lines agree, leading us to only show one Fe point per star. Ti I and Ti II abundances are shown 
with separate symbols. For the purposes of this discussion, we will ignore the BHB stars, as they 
present special challenges. See 114.41 for more details. Since our targets are all confirmed members 
of M68, the Teff-[Fe/H] trends displayed in the upper panel of Figure 0] (SPEC parameters) are 
unphysical and unexpected. Metallicities obtained from the adopted COMB (photometric, bottom 
panel) Teg still show a positive trend with increasing temperatures, but it is much less pronounced. 
Additionally, the COMB [Fe/H] values agree better with current metallicity estimates of M68. We 
attempted to eliminate the Teg-[Fe/H] trend in the RHB by using microturbulent velocities to 
correct for any differences between our stars. This approach, however, led to values ranging from 
3.0 km s“^ to 16.0 km s“^, far too high for RHB stars. 

Perhaps the more important reason for adopting our final approach can be seen in Figure O 
where we plot Alog( 5 ) {= logig)pHo^ - log(fi')comp„ri^.on value) versus Tgg for the RGB and RHB 
evolutionary stages. The top panels of each column compare spectroscopic and photometric log(( 7 ) 
values from previous studies, while the middle and bottom panels compare our SPEC and COMB 
log( 5 ) values to the PHOT ones (see M.ll for more details about photometric log( 5 ') values). Clearly, 
the differences between the purely photometric and our COMB log( 5 ') values, which are displayed 
on the y-axis in Figure O are smaller than those of any other approach and exhibit little to no 
trend with increasing Tgg. These results will be discussed in more detail in 114.31 

As a final note, we call attention to stars 117, 160 (RGBs) and 324 (BHB), for which we were 
unable to derive atmospheric parameters using either the SPEC or COMB approaches outlined 
above. Lee05, who also analyzed stars 117 and 160, believe that star 117 is an AGB rather than 
a RGB star. However, we suggest that both 117 and 160 might be extreme examples of RGB tip 
stars, which are inherently difficult to analyze using standard spectroscopic methods. For these 
two stars, the adopted Tgg and log(( 7 ) were obtained using the the purely photometric approach 
described in 114.11 For and [Fe/H], the average of all other RGB stars was chosen. Despite 


these differences in assigning atmospheric parameters, the resulting abundance patterns of these 
stars are in good agreement with other stars in our sample. If these stars are, in fact, AGB stars, 
our assumptions about the their atmospheric structures could be false, explaining our analytical 
difficulties. Star 324 (BHB) exhibited an unusually low S/N, which, in combination with it being 
a BHB star, resulted in less than 20 measurable Fe lines, making a spectroscopic analysis very 
difficult. For this star, we performed 4 iterations of the COMB approach and then adopted the 
resulting parameters. 


4.3. Further Motivations for Our Atmospheric Parameter Derivation Approach 


In the previous section, we highlighted two reasons for adopting our atmospheric parameter 
derivation methodology. In addition to quantifying these advantages in this section, we will also 
contrast our COMB parameters to literature values in the following order: Teg, log( 5 ), and 
finally [Fe/H]. 


Figure [6] compares the COMB and SPEC Teg values listed in Table 0] for the RCBs and RHBs. 
Teff values we derived from photometry are systematically higher than spectroscopic ones: (ATeg) 
= {TeSPHOT ~TeSSPEc) = 159 (u = 80) K for our study. Lee05 found a similar systematic upward 
shift of photometric temperatures, with their data giving (ATeg) = 100 {a = 61) K. For the two 
stars shared by the two studies (117 and 160), differences between their photometric and our COMB 
Tegvalue are 60 K (star 117) and 32 K (star 160). These small differences a re most likely caused 
by LeeOS’s usage of slightly different IRFM calibrations (jAlonso et al.lll999l ). as well as different 
(m — M)v and E{B — V) values. Overall, the results of the two studies are comparable. 


For RHB stars in M68, there exists no previous high-resolution spectroscopic literature refer¬ 
ence. We will therefore compare our results to the M15 RHB results of Sobll. The average offset 
between the photometric and spectroscopic Teg values of our study is (ATeg) = 321 (a = 146) K. 
Sobll found a much smaller average difference of (ATeg) = 51 (a = 272) K. The constant offset 
between the photometric and spectroscopic Teg values of our study (see Figure [6]) suggests that 
perhaps our adopted reddening value was a bit too high, causing hotter photometric temperatures. 
Unfortunately, there are no recent M68 RHB studies available that allow us to further explore this 
difference. 


A far greater discrepancy between photometric and spectroscopic parameters is present in 
derived log(g) values. The RCB side of Figure [5] clearly shows a trend between adopted final Teg 
and Alog(g) values (dehnition given above). In the upper panel of this hgure, we have included 
the Lee05 Alog( 5 ) values, which exhibit a slight upward trend of Alog( 5 ) with increasing (B — V) 
temperature. Our data (middle and bottom panel) shows a similar trend. In the specific case of star 
160, Lee05 derive a spectroscopic log( 5 ) of 0.0 dex, whereas we were unable to derive a spectroscopic 
log( 5 ) value. For the photometric log(g), LeeOS derive a value of 0.7 dex, close to our value of 0.65 
dex. Civen difference in adopted distance modulus and reddening value, this discrepancy is not 
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serious. For star 117, Lee05 derive a spectroscopic log(( 7 ) of 0.3 dex. We were again unable to 
derive a spectroscopic log( 5 ). The photometric log( 5 ) values are very close, with ours being 0.75 
and Lee05 giving a value of 0.8. The average difference between the photometric and spectroscopic 
log( 5 ) obtained by Lee05 is Alog( 5 ) = 0.67 (cr = 0.10), while the difference between our photometric 
and COMB log( 5 ) values is only Alog( 5 ) = 0.47 (cr = 0.18). 


The purely photometric study of ICarretta et al.l (2009, hereafter Car09) also has two RGBs in 
common with our study, 57 and 79. Since Car09 employed nearly identical {m — M)v and E{B — V) 
values, differences between our COMB and their final parameters are: -43 K in Tgfj and -0.31 in 
log( 5 ) for star 57, while star 79 exhibits differences of -29 K in Tgff and -0.44 in log(( 7 ). Given the 
difficulties of analyzing extreme RGB tip stars such as 160 and 117, we lend more weight to the 
differences between our study and Car09, which, as just demonstrated, are not severe. 


The Alog( 5 ') comparisons for our RHB stars are shown on the right hand side of Figure [5j The 
top panel shows the Sobll M15 data, which seems to exhibit a fairly strong Alog(( 7 )-Teff trend if 
compared to our data (middle and bottom panel). Our COMB Alog( 5 ') values give (Alog(g')) = 
0.22 (cr = 0.09), while (Alog(( 7 )) = 0.72 (cr = 0.41) for our SPEC parameters. The Sobll data 
exhibit (Alog( 5 )) = 0.36 (cr = 0.38). Clearly, the average differences (dashed lines) as well as 
Alog( 5 )-Teff trends are minimized for our COMB parameters in both the RGB and RHB stars. 


As allu ded to in H3.ll all of our d erived log( 5 ') values would be ~ 0.2 dex lower if we had 
adopted the iMelendez fc Barbuvl ([20091) based Fe II log{gf) values, instead of the NIST values. 
Hence, the Alog( 5 ') values in Figure [5] would be enhanced by about 0.2 dex. However, even with 
such an enhanced difference, the current Tgff trends would still exist. An additional effect that we 
have neglected so far is the temperature dependence of the physical log( 5 ) equation given in ^4.11 
To quantify this effect, we adopted our SPEC temperatures to determine photometric log( 5 ') values, 
which lowered all of our physical log(( 7 ) values by ~ 0.1 dex. Since this value is much lower than 
our typical log(( 7 ) uncertainties (see Table [5]), we can safely disregard the Tgff dependence of the 
photometric log( 5 ). In fact, Tgfj changes of 600 K or more are needed to reproduce log(g') shifts 
equivalent to our derived log(( 7 ) uncertainties. Please see H4.6l for more details. 


Figure!^ shows a co mparison of our derived to those of previous studies and the calibrations 
of iGratton et al.l (jl996l l. In particular, we used the theoretical PARSEC and log( 5 ) values, 
applied these to the RGB and RHB calibrations of iGratton et al.l (jl996l ). and plotted the results 
as thick dashed lines in Figure [3 with the l abel ‘PARSEC Our RGB results agree well with 
those of previous studies (jCavrel et al.l 120041 ) and the RGB calibrations in the high temperature 
end. At low tei nperatures (Tpg ~ 42 00 K), our microturbulent velocities seem to deviate from the 
empirical fits of iGratton et al.l (Il996l ). However, since stars in this temperature regime are difficult 
to analyze, we do not lend much weight to this difference. Individual comparisons with previous 
studies of Lee05 and Gar09 are not possible since Gar09, being a purely photometric study, do not 
give values and for the two stars shared with Lee05 (116 and 170), we adopted average Et values. 
See ^4.21 for more details. Our RHB stars agree well with those of iFor fc SnedenI (j2010l l. The 
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P ARSEC trend of decr easing with decreasing Teg is also followed by both our sample and that 
of iFor &: SnedenI (j201Cll b Given this good agreement, we regard our RHB values as satisfactory. 
We defer discussion of BHB values to the next section. 



between our and previous studies. 


All of the atmospheric parameters described above were derived without considering Fe NLTE 
effects. Our assumption of ionization balance, which we used to obtain logfq) values, has b een 


identified as potentially leading to incorrect atmospheric parameters (i.e. iFabrizio et al.ll2012l and 


Bergemann et al.M2012ll. To quanti fy the severity of these effects in our sample stars, we used 


Figure 4 of iBergemann et al.l (|2012l i to estimate NLTE corrections for our derived log(g) values. 
Unfortunately, we were only able to do this for some of our RGBs (all except 117, 160, 450, 
472, and 481). We are currently unaware of any published NLTE calculations for RHB and BHB 
stars. After applying the NLTE corrections to our RGBs, we found that our log( 5 ) values were 
raised by approximately 0.3 dex, which, in turn, would essentially eradicate the difference between 
our PHOT and GOMB log( 5 ) values shown in Figure [5l However, it would also lead to greater 
differences between our atmospheric parameters and PARSEC isochrones. This will be discussed in 
further detail in ^4.51 Since NLTE Fe calculations are only available for a small selection of RGBs 
in our sample, we decided to ignore these effects for our analysis (see ^ for a discussion of the 
effects of this choice on derived abundances). A future effort that considers NLTE effects over the 
whole parameter range of evolved stars in this cluster is welcome. 


4.4. Challenges of the BHB stars 


The BHB stars in our sample suffer from photometric and spectroscopic deficiencies, which 
need to be discussed before proceeding to analyze their abundances. Photometric difficulties include 
the following: 


(I) the lRamfrez &: Melendez! (j2005l ) calibrations for our metallicities become unreliable at ap¬ 
proximately 7000K due to a lack of reliable calibration stars 


(H) at Teff > 7000K, the peak of the stellar blackbody curve lies at shorter wavelengths than 
the centers of the B ( 4400 A) and V ( 5500 A) bandpasses. Therefore, fluxes are now being 


— [F6/H]previous study “ [Fc/Hjour study 
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measured in the temperature-insensitive Rayleigh-Jeans tail of the Planck distribution, resulting in 
Tgff values that are extremely insensitive to changes in either the B oi V flux. Other photometric 
fluxes (mainly K) are available for some of our stars, but only (B-V) is available for the entire 
sample. Moreover, the offset between the flux peaks of the BHB stars and V-K bandpasses is 
worse than for B-V. We therefore did not consider any photometric Tgu or log(( 7 ) values derived 
from (V-K). 


The spectroscopic temperatures are also afflicted by certain weaknesses. Due to the higher 
temperatures of the BHB stars, many lines, especially those with high excitation potentials, are 
not present in their stellar spectra. The reason for this behavior is that the strength of absorption 
lines is dependent on both the line and continuous absorption. The continuous absorption, mainly 
due to H~ (and neutral H in the violet and near-UV spectral regions), increases sharply with 
temperature and is much larger in BHB stars than in RHB and RGB stars. In our BHB stars, 
we therefore see a severe drop-off in the availability of weak high EP lines. The absence of these 
lines introduces a bias when determining the atmospheric parameters purely from spectroscopic 
line analysis. Furthermore, the increased temperatures also eliminated one of our atmospheric 
parameter diagnostic elements, Ti I. The absorption lines of this species are not strong even in 
RHB stars, and become undetectable in the BHB stars. Given the difficulties in both photometric 
and spectroscopic approaches, we advise the reader to treat all of the atmospheric parameters of 
the BHBs with caution. A more involved treatment of M68 BHB stars can be found in lBehn (j2003l b 


Our simple analysis of the BHBs, however, seems to produce values that are somewhat 
comparable to those of previous studi es. Figure [3 show s a comparison o f our d ata to a linear fit 
(thin dotted l i ne) of the RR-Lyrae of For et al. ( 2011 1 and Govea et al. ( 2014 1 and the BHB of 


For &: SnedenI (1201011. For this fit, w e have excluded BHB stars with v sin i > 15 km/s. It has 


been suggested bv iGovea et al.l (|2014l l that BHB stars with larger rotational velocities suffer from 
abnormally large values; see that study for more details on this point. Our values fit very well 
with all of the previous data and we therefore regard our microturbulent velocities as satisfactory. 


4.5. Comparison of Final Atmospheric Parameters with Isochrones 


PARSEC isochrones (jBressan et al.M2012l ) exhibit several attributes which ultimately led us 
to choose them as the master isochrones for this project. Most importantly, PARSEC isochrones 
allow for the computation of isochrones for arbitrary input metallicities and ages. This is due to the 
usage of two different types of opacities during t he calculation s. For the low temperature regime, 
opacities are obtained from the ./ FSOPUS code (Marigcj booil l. while high temperature opacities 
are calculated using OPAL 1996 ( Iglesias &: Rogerslll996l l data. Additionally, the latest version of 
the freely available FREEEO^ is used to derive the equations of state. Another very important 


®littp: //freeeos.sourceforge.net / 
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aspect of the PARSEC isochrones show the - log{g) relationship beyond the tip of the RGB, 
a usual stopping point for other sets of calculations. All of these improvements over older sets of 
isochrones should provide a good estimate of the relationship between surface gravities and effective 
temperatures for different evolutionary states in M68. 

Figure [8] compares a PARSEC isochrone calculated assuming an age of 11.2 Gyr and a metal- 
licity of Z = 1.178x10“^ (both consistent with current M68 estimates) with our COMB (filled 
symbols, usual colors) and SPEC (unfilled symbols) parameters. Overall, our final data matches 
the PARSEC isochrone much better, providing additional motivation for our COMB approach. We 
also included NLTE corrected parameters in Figure [8] (black triangles). As mentioned in 114.31 these 
values show greater discrepancies with the isochrones. To close this now bigger gap between the 
NLTE parameters and PARSEC isochrones, the metallicity of our evolutionary calculations would 
have to be increased by ~ 0.4 dex, leading to an inferred [Fe/H]M68 ~ -1.80 dex, a value completely 
at odds with previous literature studies of this cluster. We also note that isochrone calculations 
are nearly age independent at M68 metallicities, and thus this difference cannot be accounted for 
by an adjustment of this input parameter. 


A comparison of the Wal94 photometry data to this isochrone in color space is shown in 
Figure El In both of th ese figures, we have also included the latest version of BaSTI isochrones 
(jPietrinferni et al.l 120041 ) with the age and metallicity of M68. In the Tefj - log(g) plane, the 
difference between the two isochrones is negligible, while in the color plane, the two isochrones give 
very different answers. The better fit of our data to the PARSEC isochrone tracks in the color 
plane provides further motivations for using these calculations. 


4.6. Uncertainties 

Instead of deriving uncertainties for all stars in our sample, we chose representative members 
of four different evolutionary stages in our sample: lower RGB (172), RGB tip (472), RHB (36) 
and BHB (337). To obtain the uncertainty, we hrst calculated the errors in the photometric 
temperature by considering the color uncertainties given in Wal94. For all of the evolutionary 
stages, this amounts to ctphot — 20K for the RGB stars and ctphot — 70K for RHB stars. Since 
errors stemming from the scarcity of available measured lines dominate in the BHBs, this part 
of the error analysis was not performed for this evolutionary stage. In addition to photometric 
errors, we also considered errors in the reddening, which amounted to cte^b-V) — 10^ 

RGBs and ue{B-V) — 45K for RHBs. The measurement uncertainty contribution to our total 
was simulated by adjusting the temperature of stars until we achieved a shift in Fe I abundances 
equal to the standard deviation of the abundances implied by the stellar Fe I lines. The errors from 
the photometry and measurements were added in quadrature to give the final Tgfj uncertainty. 
Uncertainty in log(g) was then determined by applying the upper and lower values of Tgff to our 
stars and then again requiring ionization balance. For uncertainties, we repeated the process 
that we used to determine final atmospheric parameters (see Sec. 1121), but this time applying the 
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upper and lower temperature limits. The final uncertainties determined in this manner can be 
found in Tabled We did not consider atmospheric metallicities in our error analysis, as they have 
a negligible effect on total parameter errors. 

The differences in uncertainties between the evolutionary stages is easily explained. Atmo¬ 
spheres of RGB stars are much more sensitive to Teg and log( 5 ) changes than the more evolved 
stars due to the greater number of observed lines and greater range of excitation potentials, ion¬ 
ization state and log{gf) values. As Teg and log( 5 ) increase, we see a drop off in the number of 
observed lines and therefore less sensitivity to changes in atmospheric parameters. Hence, our 
RHB have larger uncertainties than our RGBs, but smaller uncertainties than the BHBs, just as 
expected. 


5. ABUNDANCE ANALYSIS 

In this section we present the results of our abundance analysis. Unlike the discussion of the 
atmospheric parameters, we will discuss these not individually by evolutionary stage, but rather 
attempt to give an overview of post main sequence stars in M68. All of the following abundances are 
derived using the GOMB parameters given in TableHl Table[8]gives average abundances in standard 
[X/Fe] fashion for each star in our sample, while Table [9] summarizes the average abundances for 
each evolutionary stage and M68 as a whole. The first two lines in these tables give [Fe/H] for the 
corresponding star/evolutionary stage. We note here that the Fe I and Fe II of stars II7 and 160 
do not agree. The reason for this non-agreement is explained in 114.21 

For completeness, we will now briefly discuss the sensitivites of our abundances to the adopted 
atmospheric parameters, which we determined by deriving abundances using the SPEC and PHOT 
parameters given in Table 0] for stars 472, 172, 36, and 337 (stars used to derive atmospheric 
parameter uncertainties). Since we did not derive PHOT and [Fe/H] values, we adopted those of 
the COMB approach for the stars listed above. This is justified by repeating the analysis described 
in 114.21 but this time fixing log( 5 ') to the photometric value, which showed that the resulting 
PHOT and [Fe/H] values differ very little from those of the COMB approach. The derived 
absolute abundances shifted as a result of the adoption of the PHOT and SPEC parameters, but 
the metallicity, which we use to normalize our abundances, changed by approximately the same 
amount (see Table 0]). Therefore, the resulting normalized abundances differed from the COMB 
abundances by less than the abundance uncertainties discussed in H5.3I and listed in Table 01 We 
thus anticipate no serious changes to our derived abundances due to the adoption of either the 
PHOT or SPEC parameters. 

A summary of the element ratios ([X/Fe]) for 8 < Z < 70 in box-plot form can be found in 
Figures [ini[iii and 071 Unfortunately, the quality of our spectra did not allow us to make any 
inferences about differences in He abundances between MIS and M68 (see [|T]for more details). As 
before, red data represents the RGB, yellow data represents the RHB, and the blue data shows the 
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BHBs. We have also added data from previous studies, which is shown as grey boxes. In these 
box plots, the mean elemental abundance is represented by the central bar, while the upper and 
lower edges of each box give the first (lower edge) and third (upper edge) quartile of the plotted 
data. The whiskers constitute a ~ 3cr boundary. Outlier data beyond the whiskers is depicted 
as solid symbols with corresponding colors. Table [7] gives a summary of previous studies used for 
comparison in Figures m\iii and dll A list of lines, as well as the method used to derive the 
abundances (EW or synthesis) are given in Table [3l 


5.1. Light and Fe-Group Elements 


Carbon and Nitrogen: Due to the low metallicity of M68 and the S/N limitations of our spec¬ 
tra, determining C I and N I abundances for our sample stars was very difficult. After co-adding 
the spectra of our four coolest RGB stars, we were able to estimate (^^C/^^C) ~ 5—7 for the coolest 
giants in M68. The quality of our stellar spectra precluded determination of values for individual 
stars. Better quality data are needed to make a more accurate carbon isotopic assessments. How¬ 
ever, these results are in line with the low values obtaine d for other globular clusters such 

as M3 « 6. Pilachowski et al. 2003 1. M4 ^ 5, Brown fc Wallerstein 1989 ). and 

M22 « 6, Smith fc SuntzefJ 1989 ). Assuming a value of 6, we then were able 

to estimate [C/Fe] ~ -0.5, and [N/Fe] ~ 1. Figure fT2] shows the co-added RGB tip spectra and a 
corresponding synthesis for the CN region around 3880 A. Noise limitations are evident, even in 
these co-added spectra. C I and N I abundances will not be further considered in this paper given 
the difficulties in trying to determine these abundances. 

Oxygen and Sodium: The O I abundances for the RGB stars in our sample were determined 
by synthesis of the [O I] 6300.3 A line, while for three of our RHB and all of the BHB stars, 
oxygen abundances were determined from the O I triplet at 7771.94 A, 7774.17 A, and 7775.39 A. 
NLTE corrections are negligible for any forbid den transitions, but such effects can be quite large for 
abundances derived from the O I triplet (eg., ISitnova et al.ll2013l and references therein). We used 
the estimates given in their Table 11 to extrapolate NLTE corrections for our final temperatures 
by deriving a logarithmic relation between these corrections and temperature, which is given by: 


[0/Ee]NLTE — [0/Ee]LTE = —1.43 loge(refj) -I- 12.066. 

NLTE corrections for our stars range from —0.35 dex to —0.76 dex. 

Na I abundances in our stars were derived using EW measurements from four Na I lines: 
5889.95 A, 5895.92 A (the D lines), 8183.26 A, and 8194.82 A. The EWs for these lines can be 
found in Table [3l We der ived NLTE corrections using the INSPEGT0 web interface, which is 
based on iLind et al.l (120111 ). Unfortunately, the parameters of some of our stars exceed the T^s- 
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log( 5 ) limits. For such parameters, we extrapolated the NLTE corrections given by the websites to 
our specihc log(g) and Teg combinations. In the case of RGB stars, the limit was set by log(( 7 ) as 
the website did not accept values lower than 1.00. For the BHB stars, the limit was set by Teg. 
NLTE corrections obtained is this manner should not be viewed as definite, but rather as zeroth 
order estimates. 


Erom Eigure [TOl it seems that O I abundances are enhanced in our RGB stars as compared 
to other evolutionary stages and previous studies. However, since we are using different features 
to determine the O I abundances in different evolutionary groups, caution is warranted when 
interpreting this abundance trend. In Figure [T3l we compare Na I and O I abundances. The anti¬ 
correlation of t hese two elements, w hich seems to be exclusive to GGs, has been confirmed by many 
studies (Gar09. iGratton et al.ll2n04l and references therein). It is currently be lieved that this tren d 
results from pollution of the cluster medium by an earlier generation of stars (jGratton et al.ll2nn4l ). 
Unfortunately, we cannot confidently say that we observe such an anti-correlation in our sample, 
since we were only able to derive both Na I and O I abundances for 11 out of our 25 stars. However, 
the abundances derived for these stars fall within the same general ranges, perhaps slightly elevated, 
as those of Gar09 (grey symbols), which consist of some of the most metal-poor GGs known: NGG 
6397 ([Fe/H] -2.02), M55 ([Fe/H] -1.94 ). M15 ([Fe/H ] « -2.37), and M30 ([Fe/H] « -2.27). 

All quoted metallicities were obtained from iHarrisI (jl996l ). Data for stars 117 and 160 (see 114.21 
for more details) are shown in green again. These two stars seem to have particularly high [0/Fe] 
values, which, in light of the problems associated with deriving their atmospheric parameters, are 
no cause for serious concern. We conclude that, on average, M68 is not overabundant in sodium, 
but our data does suggest an oxygen overabundance. Additionally, if our data is combined with the 
Gar09, we clearly see an anti-correlation of our derived oxygen and sodium abundances. A closer 
look at Figure [l3] also appears to reveal a systematic difference of O I abundances between the 
RG B, RHB, and BHB. Such differences have been observed and analyzed by several other authors 
(see iMarino et al.l l201ll and references therein). Given that we employed different lines to obtain 
our abundances of each of the evolutionary stages and relied on crude corrections to account for 
NLTE effects, we will not attempt to dissect the details of these differences here. 


We now compare our oxygen abundance results with those of Lee05 and Gar09. Gar09 derived 
similar oxygen abundances also using the 6300.3 A [O I] transition for several different GGs. The 
average metallicity for the Gar09 M68 stars is slightly higher than ours, and since they adopted 
purely photometric log( 5 ') and Tgg values, their slightly lower O I abundances are expected. In the 
specihc case of the two overlapping RGBs, the differences are A[O/Ee]570 = —0.35 for star 57, 
while for star 79 we could not derive an O I abundance. The discrepancy between our results and 
those of Lee05 are more difficult to understand. Eor star 160, A[0/Fe]i6o = —0.43 and for star 117 
A[0/Feii7] = —0.53. However, Lee05 employed a very different atmospheric parameter derivation 
methodology. Eurthermore, they did not force ionization balance between neutral and ionized Ee 
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lines, resulting in Fe II abundances that are on average, 0.37 dex higher than the corresponding Fe I 
abundances. They chose to normalize their oxygen abundances using these elevated Fe II values, 
naturally leading to higher [0/Fe] ratios. 

The differences between our sodium abundances and those of Car09 are A[Na/Fe ]57 = +0.17, 
and A[Na/Fe ]79 = +0.14 . These differences can mostly be attributed to the slightly different param¬ 
eter derivation methodology and NLTE correction algorithm employed by CarOO. The discrepancies 
between our work and that of LeeOS are a bit greater: A[Na/Fe]ii 7 = —0.27 and A[Na/Fe]i6o = 
+0.36. Again, these differences, considering the disparate approaches, are negligible. We note here 
that Na I NLTE corrections could not be derived for star 160 using the methods described above. 


Aluminum: All of our aluminum abundances were derived from the A1 I resonance lines 

at 3944.01 A and 3961.53 A. Due to severe NLTE effects, these lines are known to give abun- 
dances which are, on average, a factor of 6 lower than those derived from other A1 I lines (e.g., 
Baumiiller &: Gehrenlll997l . lAndrievskv et al.ll2008lL Unfortunate ly, these are the only measurable 
A11 lines in our spectra. In their Figure 2. lAndrievskv et al.ll2008l give NLTE corrections applicable 
to some of our giants. We estimate corrections of 0.4 dex for stars 79 and 172, a correction of 0.55 
dex for stars 226 and 256, and a correction of 0.5 for star 44 0. Atmospheric pa r amete rs for all other 
stars in our sample are outside the recommended limits of I Andrievsky et al.l (j2008l ) and therefore 
are not considered for NLTE corrections. Table [6] gives the NLTE A1 I abundances, as well as 
NLTE corrected Na I abundances. Clearly, there is a spread amongst A1 I abundances in M68. 
However, all low/high [Al/Fe] ratios correspond to low/high [Na/Ee] ratios. Figure fT^ shows this 


well- know n correlation b etween [Al/Fe] and [Na/Fe] for globular clusters (e.g Car09, llvans et al 


2001 


and Shetrone 1996h . In the interest of having the maximum number of data points available, 


we decided to use the Al I LTE abundances. Since the NLTE corrections are almost constant for 
our stars (see discussion above and Table [6|), there is no danger in destroying any abundance corre¬ 
lations by using the LTE abundances. The clearly positive correlation between these two elements 
provides evidence for primordial abundance enhancements by earlier generations of stars. We will 
not atten ipt to explore these primordial abundance enhancements using our data, but refer the 
reader to lGratton et al.l (|200J) for more information. 


Lee05, who derive ([Al/Fe]) = 1.08, and Car09, who derive ([Al/Fe]) < 0.74, measured their 
abundances using the subordinate Al I lines at 6696.03 A and 6698.67 A. We note that Car09 
explicitly sta te that all their A l I ab undances are upper limits, which, in combination with the 
statement bv lAndrievskv et al.l (j2008l l that these subordinate lines should not be visible at [Fe/H] 
< —2.5, leads us to lend less significance to the differences between our study and Lee05. We derive 
([Al/Fe]NLTE) = 0.31, which given the difficulties explained above, seems in concord with CarOO. 
To ensure that our lack of detection of these subordinate lines is not caused by noise limitations, we 
co-added RGB spectra of the stars with the highest Al I abundances implied by the resonance lines. 
Despite these efforts, we were not able to measure the subordinate lines, or Al I line at 7836.14 A, 
which suggests that [Al/Fe] < 0.78. 
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a Elements: The a and a-like elements considered in this study are Mg I, Si I, Ca I, and 
Ti l/Ti II. Their abundances, which were obtained from EW measurements, are compared in 
Figure [TOl We will discuss them in order of increasing Z. 


Mg I abundances for all evolutionary stages show good agreement, which, in combination with 
Figure fT^ where we see the expected [Mg/Fe] —[Al/Fe] anti-correlation, supplies further evidence for 
primordial abundance variations in GC stars. Dissecting the details of these abundance behaviors 
is beyond the scope of this paper. Here, we mainly demonstrate that our atmospheric parameter 
derivation approach reproduces all the classical abundance hallmarks of any GC study. Star 202, 
with [A1 l/Fe] ~ 0.4 in Figure [15] is the hottest RHB in our sample, and with Tgfj = 6257, very 
close to the RR-Lyrae gap. Therefore, the A11 abundances for this star could be compromised. We 
have also chosen not to include BHB Mg abundances in Figure [15] due to reasons listed in 1 14.41 


Car09, Lee05 and ]Behi] (|2003] i also derived Mg I abundances for their stars. They differ from 
our measurements as follows: ^mqIFc, 57 = -0.044, = -0.083, i^MgjFe ,\\7 = -0.25, and 

‘^Mg/Fe,i 60 = -0.01. All of these ratios show excellent agreement, except those of star 117. The 
large discrepancy for this star, however, is no real surprise given the hugely different approaches 
for determining atmospheric parameters and the fact that this particular star might be an example 
of an extreme RGB tip star, making its analy sis very diffi cult. We will not compare the individual 
stellar abundances for star 324 with those of iBe. 
atmospheric abundances for this particular star. 


Ill (]2003]i. s ince we were unable to derive proper 


Behr] (j2003l ) derive a ([Mg/Fe]) = 0.18, while we 


obtain ([ Mg/Fe|) = 0 .35 for our BHBs. This difference can be explained by the different approach 
taken by Behr ( 20031 ) to derive their atmospheric parameters. 


Unlike the [Mg/Fe] ratios, our derived [Si/Fe] abundances show large discrepancies between 
the evolutionary stages. This is be cause, as several stu dies have shown, Si I abundances drop with 
increasing Tefj (e.g. Figure 10 of iPreston et al.ll2006l and references therein). In the interest of 
consistency, we derived all of our Si I abundances from the 3905.53 A line, which is saturated in 
the RGB, leading to possibly unreliable abundances. If we had used the 4102.94 A Si I line for this 
evolutionary stage, our RGB abundances would be elevated by about 0.4 dex. Unfortunately, the 
4102.94 A line is not measurable in the RHBs or BHBs. 


For the RGB branch, our study yields ([Si l/Fe]RGB) = 0.40 (see Table [^, which, if com¬ 
pared to ([Si I/Fe]R,GB) values fro m other metal poo r clusters: ([Si I/Fe]RGB,Mi 5 ) = 0.55, and 


([Si_l/Fe]RGB,M92) = 0.59 (both from Sneden et al.lboOC i. seems slightly too low. However, 


Sneden et al. 


2000 use only the 5948.55 A Si I line, and a different set of gf values. Similar to ]Sneden et al 


([20001), Lee05 also derive an analogously high Si I abundance, ([Si l/Fe]RGB,M68)= 0.71, which they 
base on EW measurements of 6-7 mA for all their stars. Given that their S/N is comparable to 
ours, noise limitations could contribute to their derived Si overabundances. To ensure that our 
EW measurements are not at fault, we calculated the EW values needed to reproduce the Lee05 
abundance for several different high excitation Si I lines, including those used by Lee05. Subse¬ 
quently, we inspected the spectra of our stars for these lines and tried to measure the EW needed 
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to reproduce the Lee05 abundances. We were unable to measure or visually locate any of these 
lines in our spectra, confirming that M68 most likely does not exhibit an Si overabundance. In 
their study, Car09 derive ([Si/Fe]) = 0.40, in good agreement with our derived value. 


Our RHB Si measurements yield ([Si/Fe]) = 0.16, which is in accord with iPreston et al.l ([200a), 
who derive ([Si l/Fe]RHB,Mi 5 ) = 0.32. While we expect lower abundances with increasing (see 
above), a drop of 0.5 dex in [Si/Fe] between our average RHB abundances and those of Lee05 would 
be outside of the usual range found by previous studies. 


As a final remark on Si I abundances, we perform a star-by-star comparison with CarOO, 
which yields these differences: Asi/Fe ,57 = —0.23, Asi/Pejg = —0.051, The reasons for the large 
differences are explained above and can also be attributed to the different methodologies employed 
for deriving these abundances and associated atmospheric parameters. We were unable to derive 
Si abundances for either stars 117 or 160, precluding us from a comparison with Lee05. 


Perhaps a more important and reliable diagnostic than Si is Ca, which exhibits more measurable 
transitions and has been shown to give constant abundances for GCs with [Fe/H] < -1.00. Our 
an alysis yields ([Ca/Fe ]Mfis) = 0.36, very much in line with previous studies shown in Figure 4 
of Grafton et al. ( 2004l h They derive ([Ca/Fe]) =4-0.25 {a = 0.02) for 28 clusters with [Fe/H] < 
-1.00. A comparison with Lee05, who also derive Gal abundances yields the following: Aca/Fe,n 7 
= -0.16, and im = -0.04. For a thorough review on Ga abundances in GGs, please see 

Grafton et al. (2004). 


As a final a element, we measured the abundances of Ti I and Ti II, whose abundances provide 
a good diagnostic of possible over-ionization in stellar atmospheres. We will not compare our Ti 
abundances to those of Lee05, since the methodology of the two studies are too disparate. A quick 
look at Figure dO] reveals over-ionization in RGB stars (Ti II abundances are, on average, 0.34 
dex higher), while the agreement between Ti I and Ti H is excellent for RHB stars. For reasons 
explained above, there are no measurable Ti I lines in our BHB stars. Overall, our data shows an 
a elemen t enh a ncem ent, as expected from globular clusters. A similar enhancement in M68 was 
found by Behr ( 2003 1. 


Fe-peak Elements: Sc II, Cr I, Gr II, Mn I, Mn II, and Go I make up the Fe-peak elements 
measured in this study. Their combined mean abundances are ([X/Fe]i?e peak) = O.II (a = 0.30). 
Figures fTOl and fTTl show the behavior for these elements amongst different evolutionary stages. Sc II 
abundances agree for all our evolutionary states and the Lee05 data. The large difference between 
Cr I and Cr II abundances could be a cause for concern, since we required ionization equilibrium 
for our abundances. However, such behavior has been observed in other cluster as well as field 


stars and is currently unexplained (see iPreston et al.ll2006l and references therein). In Figure [TTl 
we show a comparison with the M92 data of iLanger et al.l (|l998l l. Clearly, their average abundance 
is much lower than ours. Due to the known problems with this element, however, we will not lend 
too much weight to this difference. 


Mn I and Mn II abundances also seem to show over-ionization and large variations amongst 
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individual stars. However, nearly all of these abundances are based on a single stellar line mea¬ 
surement. We therefore urge the reader to treat these abundances with caution. Co I abundances 
show agreem e nt be tween the RGB and RHB branches, but are again much lower than those of 


hanger et al.l (|l998l i. The reasons for this behavior are unknown. 


Lee05 measured Mn I and Sc II, and for these species our differences are: ^Sc/Fe,in = -0.24, 
-^Sc/Fe ,160 = -0.26, ^Mn/Fe,in = 0.03, and ^Mn/Fe,im = -0.04. Unfortunately, no other abundance 
comparisons are possible. 

Copper and Zinc: The last two elements in Figure [TT] are Cu I and Zn I. Lee05 also derived 
Cu I abundances and the difference between the overlapping star 160 is ^cu/Fe,im = -0.55. These 
differences can be attributed to the disparate approaches in deriving atmospheric parameters. For 
Zn I, no previous references exist. 


5.2. Heavy Elements 


n-capture Elements: Figure [T7] shows box plots of the distribution of n-capture elements in 
our stars. Unlike many previous studies, the violet extent of our spectral coverage allows us to 
include rare-earth-element species such as Dy II and Yb II. The distributio n of n-capture elem ents 

As 


can be used to infer r- or s-process enrichment (see Figures 1, 2, and 10 in ISneden et al.ll200^ 
part of r-process enrichment, we should observe low [Ba/Eu] and [La/Eu] ratios, both of which are 
clear ly present in our s ample. We therefore conclude that M68, much like other metal poor GCs 
(e.g. ISobeck et al.ll201ll ). exhibits r-process enrichment. The constant offset of ~ 0.26 dex between 
the RGB and RHB evolutionary stages can be attributed to the usage of different MOOG versions 
(see 114.21 for more details), as well as general difficulties associated with the size of the parameter 
space covered by this study. 


Figure [17] also demonstrates discrepancies between Sr II and Ba II abundances amongst the 
evolutionary stages. These are easily explained. For Sr II, the employed lines at 4077.71 A and 
4215.52 A are saturated in almost all RGB and RHB stars, resulting in unreliable abundances. 
In the hottest RHB and all BHB stars, these Sr II lines are unsaturated and yield consistent 
abundances. The Ba II discrepancies can be traced back to the usage of different sets of lines for 
obtaining abundances. The well known 4554.03 A Ba II line is saturated in our RGB and RHB 
stars, and we used the 5853.69 A, 6141.73 A, 6496.91 A in these stars. In our BHB stars, these 
Ba II lines were not available, so we reverted to the 4554.03 A line. 


Gar09 did not derive any n-capture process abundances, and therefore we will compare all of 
our abundances to those of Lee05. The differences for star 117 are /S.Ba ii/Fe,in = -0.28 (only 
element available for comparison), while those for star 160 are ///Fe,i60 = -0.29, ///Fe,i60 

= -0.25, and Aeu ii/Fe,i 60 = -0-14. Like before, A^/Fe is defined to be [X/Fe\previous study - 
[A/Fe]this study Possible reasons for the differences between our abundances and those derived by 
Lee05 have been discussed in N5.1I 
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5.3. Abundance Uncertainties 

Much like atmospheric parameter uncertainties, we derived representative abundance a values 
for each evolutionary state in our sample. The representative stars remain the same: 472 (RGB tip), 
172 (lower RGB), 36 (RHB), and 337 (BHB). We note that we derived Cr I abundance uncertainties 
for the BHBs using star 289, since we could not measure any Gr I lines in star 337. 

To obtain uncertainties for elements whose abundances are based on EW measurements, we 
took the following approach; Uncertainties for any element for which more than 3 lines were 
measured, the standard deviation of the abundances implied by the measured lines was used. 
For any element with less than 3 lines, we re-measured the EW now considering factors such 
as continuum placement and smoothing. Using these re-measured EW values, we derived new 
abundances, which allowed us to calculate uncertainties for the corresponding species. This method 
led to uncertainties of ~ 0.13 dex for most elements. 

For uncertainties for elements where spectral synthesis was used to obtain abundances, we 
re-synthesized the spectra and again considered factors such as continuum placement, smoothing 
and assumed abundance. Just like with our EW uncertainties, this method leads to final a values 
of ~ 0.13 dex. The final uncertainties obtained in the above described fashions can be found in 
Table [TOl It is also mentioned here that our uncertainties in Tgff, log( 5 '), and lead to additional 
uncertainties of 0.2 - 0.3 dex, which are not included in the results of Table fTOl 


6. CONCLUSIONS 

In this paper, we have explored the atmospheric parameters and detailed chemical compositions 
of 25 evolved members of M68. Particular attention has been paid to comparisons of the assets and 
liabilities of photometrically based and spectroscopically based parameters. From the discussion 
in ^4.21 and the evidence presented in 114.31 it is now clear that difficulties in deriving atmospheric 
parameters in an internally consistent manner over a parameter space that covers ~ 4000 K in Tgff, 
and ~ 3.00 dex in log( 5 ) exist and need to be accounted for. Luckily, for M68, as for almost any 
cluster, we have reliable reddening and distance moduli that allowed us to treat these problems by 
creating a hybrid spectroscopy-photometry approach. 

Some of the discovered weaknesses of the two standard approaches include photometric deriva¬ 
tion of log( 5 ') values. As mentioned in Sec. 14.11 the physical log( 5 ) depends on the mass of the star 
in consideration. Unfortunately, precise mass loss between the turn-off, RGB and later evolutionary 
stages is still uncertain. Our method of applying a single mass to all stars in our sample provides a 
good zeroth order approximation of the physical log( 5 ), but for increasingly detailed studies, more 
accurate stellar masses should be adopted. Spectroscopic weaknesses include scarcity of lines in 
hotter RHB and BHB stars, as well as serious metallicity trends (see Figure |4]). A combination of 
both the photometric and spectroscopic methods allowed us to address most of these problems. One 
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of the main highlights of our final adopted approach is the lack of any obvious metallicity difference 
between evolutionary stages (see Figured] or ([Fe/H]) values for RGBs and RHBs in Table [9]). Our 
analysis results in u^pe/H]) = 0.14, for all stars, including the BHBs. We remind the reader that this 
lack of metallicity difference is not inherent to our analysis or our atmospheric parameter derivation 
method. Instead, it naturally grows out of our adoption of photometric temperatures, making this 
result even more remarkable. 

In addition to the metallicities being in good agreement. Table [9] shows that the abundances 
for all other elements agree across the different evolutionary stages. For the elements that seem to 
exhibit any discrepancies, a valid and detailed explanation is given in We can therefore conclude 
that even though we adopted a non-standard, hybrid approach to deriving our atmospheric param¬ 
eters, the resulting abundances are what one would expect from a classical purely spectroscopic 
analysis. Moreover, we were also able to reproduce all classical hallmarks of GC populations, such 
as the [Na/Fe] —[0/Fe] and [Mg/Fe] —[Al/Fe] anti-correlations, Ga and Si abundances that agree 
with previous studies, as well as slight r-process enrichment. In the absence of very detailed NLTE 
calculations and/or 3D model atmospheres, such a hybrid approach may be necessary for us to 
further develop our understanding of cluster stars, at least from a stellar atmospheric perspective. 


We thank Sloan Simmons for initial work on this project, as well as Karin Lind and Jeremy King 
for helpful discussions. We also thank our referee for providing valuable suggestions and thereby 
improving the quality of this paper. Financial support for this research from the US National 
Science Foundation (grant AST-1211585) and the Rex G. Baker, Jr. Centennial Endowment to the 
University of Texas are gratefully acknowledged. 
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Wavelength [A] 


Fig. 1. — Sample spectra for all evolutionary stages in our sample. In addition to the name and evolutionary state 
of the star, we also list and log(g). This specihc wavelength region covers the CH G-band. Measuring this band 
in the RHB and BHB is nearly impossible due to their higher Tefj and log(g) values. 
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{B-V) 


Fig. 2. — CMD of M68 using IWalkeil (1994) data. The markers and colors defined here for the RGB, RHB, and 
BHB will be used throughout the remaining plots in this paper. The green points represent non-standard stars. 
Please see 44.51 for more details. 
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Fig. 3. — Here we show the implied Tefi values for each of the color indices. For the purposes of this figure, 
ATeff = TeS{v-i) ~ Tefi{B-v)- before, RGBs are represented as triangles, RHBs as circles, and BHBs as squares. 



[X/H] 
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Fig. 4. — Comparison of ionization and metallicity effects using different atmospheric parameters. The solid black 
line represents a linear fit to [Fe/H] of RGB and RHB stars as a function of temperature. 
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T,„/1000 [K] 


Fig. 5.— Comparisons of Alog( 5 ) (= \og{g)p^Q,^ - vaiu^) values as a function of temperature for 

different studies and methodologies. The solid lines represent a linear fit to the data, while dashed lines give the 
average Alog(g) offset. 
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Fig. 6.— Comparison of photometric and spectroscopic temperatures. Our data agrees well with Lee05, but differs 
quite significantly from Sobll. 
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Te,f [K] 


Fig. 7.- 


Microturbulent velo cities as a fun c tion o f evolutionary stage. 

(Il996ll . 


studies and the calibrations of iGratton et al 


Our results agree well with previous 
In this Figure, we have used a diamond to represent the 
RR-Lyrae evolutionary stage. Solid symbols represent this study, while unfilled symbols represent comparison values 
from previous studies. See ii4.3l for more details. 
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Teff [K] 


Fig. 8.— Comparison of isochrone tracks with derived atmospheric parameters. Green markers indicate derivation 
of atmospheric parameters by photometric data only. Markers without fill indicate purely spectroscopic parameters. 
See tl4.2l for more information. 
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Fig. 9.— Reproduction of Figure [8] but in color space. The fit of the PARSEC isochrones is somewhat better 
than that of the BaSTI tracks. 







[X/Fe] 
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Fig. 10.— Comparisons of derived abundances for 8 < Z < 22. Any abundance differences between the evolutionary 
stages are explained in i)5.1l Solid markers outside the boxes represent outlier abundances. Grey boxes represent 
data from previous studies. Elements marked with have NLTE corrections applied to them. 
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Fig. 11.— Continuation of Figure fTIil but with 24 < Z < 30. Abundance differences are discussed in ii5.1l 




















Intensity 
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Fig. 12. — Synthesis of the CN region around 3880 A. Noise limitations are evident from this plot. The spectrum 
shown here is a co-addition of stars 117, 160, 472, and 481, all of which are RGB tip members. 
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Fig. 13. — Correlation between [Na/Fe] and [O/Fe]. Our data is compared to that of Car09, who derived these 
abundances for RGBs in M68, NGC 6397, M55, M15, and M30. NLTE corrections have been applied to for both O 
and Na abundances. 
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Fig. 14. — The classical comparison of [Na/Fe] to [Al/Fe], which produces the standard result for globular cluster 
studies. In this plot, NLTE corrections were applied only to Na. 







[Al/Fe] 
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Fig. 15.— Anti-correlation of [Al/Fe] and [Mg/Fe] for stars in our sample. 







[Si/Fe] 
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[Mg/Fe] 


Fig. 16. — Correlation between Si and Mg abundances. The drop-off in the Si abundances with increasing 
temperatures is well known. 







[X/Fe] 
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Fig. 17. — Abundances of n-capture elements amongst the different evolutionary stages. A slight r-process 
enrichment can be observed. Further details are given in H5.2I 
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Table 1. 


Fundamental parameters of NGC 4590 



Right ascension^ 

12 36 27.98 

Declination^ 

-22 44 38.6 

Distance^ 

10.3 kpc 

Age^ 

11.2 Gyr 

Metallicity'^ 

-2.23 

E{B-Vf 

0.05 

{m — M)v'*‘ 

15.210 


’ Goldsburv et al. I 201oll 
Sollima et alj ( 2008h 


Harris! ( 1996h . (2010 edition) 
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Table 2. Photometric parameters 


Star 

V 

{B-VJ 

{V-IJ 

RGB photometry 

57 

14.213 

0.937 

1.116 

79 

13.830 

0.993 

1.161 

117 

12.693 

1.191 

1.320 

160 

12.694 

1.319 

1.418 

172 

14.488 

0.894 

1.084 

226 

15.514 

0.781 

0.981 

256 

15.202 

0.814 

1.004 

440 

14.814 

0.842 

1.030 

450 

13.450 

1.023 

1.182 

472 

13.211 

1.103 

1.244 

481 

12.867 

1.245 

1.359 

RHB photometry 

36 

15.588 

0.524 

0.723 

47 

15.545 

0.567 

0.775 

202 

15.591 

0.446 

0.676 

334 

15.586 

0.558 

0.727 

403 

15.465 

0.586 

0.798 

454 

15.447 

0.507 

0.740 

458 

15.638 

0.503 

0.700 

533 

15.574 

0.549 

0.746 

547 

15.613 

0.459 

0.685 

BHB photometry 

170 

15.714 

0.193 

0.323 

289 

15.603 

0.264 

0.295 

324 

15.772 

0.204 

0.287 

377 

15.701 

0.212 

0.348 

391 

15.704 

0.250 

0.311 
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Table 3. EW lists 


Ion 

A (A) 

E.P. (eV) 

log(g/) 

57 

EW 

log e 

79 

EW 

log e 

117 

EW log e 

Method 

Ref. 

O I 

6300.30 

0.000 

-9.820 


7.07 


7.04 


7.15 

synth 

1 


7771.94 

9.150 

0.370 







EW 

2 


7774.17 

9.150 

0.220 







EW 

2 


7775.39 

9.150 

0.000 







EW 

2 

Na I 

5889.95 

0.000 

-0.190 



221.6 

3.97 



EW 

3 


5895.92 

0.000 

0.110 

205.5 

4.19 

194.2 

3.94 



EW 

3 


8183.26 

2.100 

0.240 

67.1 

4.21 

55.5 

4.01 

75.1 

4.06 

EW 

3 


8194.82 

2.100 

0.490 

111.9 

4.53 

82.5 

4.10 

112.6 

4.27 

EW 

3 

Mg I 

3829.36 

2.710 

-0.227 



232.6 

5.47 



EW 

4 


3832.31 

2.710 

0.125 







EW 

4 


4571.10 

0.000 

-5.623 

105.9 

5.54 

106.3 

5.41 

158.9 

5.59 

EW 

4 


5172.70 

2.710 

-0.393 





135.6 

5.86 

EW 

4 


5183.62 

2.720 

-0.167 







EW 

4 


5528.42 

4.350 

-0.498 

100.3 

5.54 

117.1 

5.74 



EW 

4 

A1 I 

3944.01 

0.010 

-0.638 



215.2 

4.16 



EW 

3 


3961.53 

0.010 

-0.336 



177.2 

3.43 



EW 

3 

Si I 

3905.53 

1.910 

-1.041 

293.6 

5.55 

255.2 

5.36 



EW 

5 

Ca I 

4226.74 

0.000 

0.244 







EW 

6 


5588.76 

2.530 

0.210 

69.1 

4.12 

73.7 

4.13 

95.1 

4.26 

EW 

6 


5857.46 

2.930 

0.230 

51.1 

4.30 

43.0 

4.14 

66.6 

4.31 

EW 

6 


Note. — The machine-readable version of the entire table is available in the online journal. A portion of 
the table is shown here for guidance concerning form and content. The abundances listed in this table are 
given for corresponding individual transitions. 
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Table 4. Atmospheric parameters 




PHOT 

SPEC 


COMB 


Star 

S/N 

Teff 

log(5) 


log(5) 

Teff 

log(5) 

[Fe/H] 

6 

RGB parameters 

57 

140 

4626 

1.54 

4468 

0.83 

4626 

1.25 

-2.51 

2.31 

79 

130 

4553 

1.35 

4418 

0.50 

4553 

0.91 

-2.52 

2.41 

117+ 

200 

4310 

0.75 



4310 

0.75 

-2.52 

2.33 

160+ 

230 

4168 

0.65 



4168 

0.65 

-2.52 

2.33 

172 

100 

4687 

1.69 

4515 

0.82 

4687 

1.26 

-2.45 

2.30 

226 

90 

4844 

2.18 

4655 

1.22 

4844 

1.71 

-2.51 

2.05 

256 

90 

4789 

2.02 

4532 

1.00 

4789 

1.67 

-2.58 

1.97 

440 

100 

4756 

1.86 

4505 

0.98 

4756 

1.69 

-2.47 

2.04 

450 

150 

4513 

1.18 

4342 

0.02 

4513 

0.51 

-2.65 

3.08 

472 

150 

4414 

1.02 

4305 

0.25 

4414 

0.69 

-2.53 

2.38 

481 

180 

4247 

0.77 

4257 

0.23 

4247 

0.18 

-2.47 

2.42 

RHB parameters 

36 

90 

5905 

2.62 

5532 

1.72 

5905 

2.42 

-2.46 

3.01 

47 

100 

5724 

2.55 

5392 

1.70 

5724 

2.35 

-2.47 

2.57 

202 

90 

6257 

2.73 

6180 

2.30 

6257 

2.46 

-2.31 

3.45 

334 

30 

5756 

2.57 

5220 

1.37 

5756 

2.40 

-2.42 

2.93 

403 

100 

5649 

2.49 

5198 

1.30 

5649 

2.20 

-2.49 

2.71 

454 

70 

5972 

2.59 

5860 

2.27 

5972 

2.53 

-2.36 

3.24 

458 

80 

6001 

2.67 

5670 

2.71 

6001 

2.30 

-2.48 

3.36 

533 

60 

5792 

2.58 

5430 

1.63 

5792 

2.33 

-2.40 

2.61 

547 

100 

6197 

2.72 

5885 

2.03 

6197 

2.58 

-2.37 

3.81 

BHB parameters 

170 

80 

7848 

3.24 

7810 

3.43 

7848 

3.50 

-2.06 

2.30 

289 

80 

7410 

3.10 

7680 

3.05 

7410 

2.59 

-2.58 

1.97 

324+ 

50 

7792 

3.25 



7792 

3.18 

-2.25 

2.17 

337 

80 

7700 

3.20 

8010 

3.68 

7700 

3.24 

-2.08 

2.02 

391 

100 

7444 

3.15 

7497 

3.11 

7444 

3.04 

-2.29 

2.39 


^These stars are represented with green markers in Figs. [2 [3 and [3 See 14.21 
for details. 
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Table 5. Parameter uncertainties 


Evolutionary stage 


'^log(g)[c“/®"] 

[km/s] 

early RGB 

100 

0.30 

0.20 

late RGB 

75 

0.34 

0.25 

RHB 

150 

0.30 

0.40 

BHB 

200 

0.35 

0.20 
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Table 6. [Al/Fe] and [Na/Fe] NLTE abundances 


Star 

[Al/Fe]NLTB 

[Na/Fe]NLTB 

79 

0.26 

0.08 

172 

-0.34 

-0.26 

226 

0.73 

0.63 

256 

0.50 

0.30 

440 

0.42 

-0.14 
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Table 7. Comparison study references 
Element Reference 


O I ^arretta_et_^ (2009) 

Na I ^arretta_et_^ (2009) 

Mg I ^arretta_et_^ (200^ 

All 

Si I ^axretta^_et_aL (2009) 

Ca I Jjee_^t_al (2005) & Jjan^er_et_aL (1998) 

Sc II ^jge^e^^ (2005) & Ijauge^_e^^ (1^98) 

Ti I 2jee_^^_al (2005) & Jjauger_et_^ (1998) 

Ti II 2jee_^^_al (2005) & Jjau^er_et_^ (1998) 

Cr I Jjangei_et_^ (1998) 

Cr II 

Mn I Jjee_et_^ (2005) 

Mn II 
Fe I 
Fe II 

Co I 2janger_et_^ (1^98) 

Cu I Jjee_et_^ (2005) 

Zn I 

Sr II Jjau^ei_e£_^ (1998) 

Y II 
Zr II 

Ba II Jjee_et_al (200^ & 2jan^er_et_^ (1998) 
La II (2005) 

Nd II 

Eu II (2005) 

Dy II 
Er II 
Yb II 
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Table 8. Average abundances for individual stars 


Ion 

57 

{[X/Fe]> 

79 

([X/Fe]> 

117 

([X/Fe]> 

160 

([X/Fe]> 

172 

([X/Fe]> 

226 

([X/Fe]> 

256 

([X/Fe]) 

Fe 1+ 

-2.51 

-2.52 

-2.39 

-2.49 

-2.45 

-2.51 

-2.58 

Fe 11+ 

-2.51 

-2.52 

-2.16 

-2.18 

-2.45 

-2.51 

-2.58 

Ol 

0.80 


0.89 

0.76 




Nal 

0.31 

0.08 

0.27 

0.10 

-0.26 

0.63 

0.30 

Mg I 

0.45 

0.46 

0.64 

0.51 

0.33 

0.00 

0.34 

All 


-0.14 



-0.74 

0.18 

-0.05 

Si I 

0.55 

0.37 



0.53 

0.47 

0.53 

Cal 

0.34 

0.33 

0.43 

0.33 

0.41 

0.36 

0.39 

Ti I 

0.09 

0.15 

0.08 

-0.08 

0.14 

0.09 

0.16 

Ti II 

0.33 

0.28 

0.72 

0.75 

0.25 

0.30 

0.41 

Sc II 

0.16 

0.12 

0.43 

0.42 

0.18 

0.10 

0.20 

Cr I 

-0.15 

-0.08 


-0.06 

-0.02 

-0.17 

0.01 

Cr II 

0.36 

0.26 

0.56 

0.54 

0.29 

0.41 

0.34 

Mnl 

-0.39 

-0.37 

-0.25 

-0.26 

-0.36 

-0.36 

-0.38 

Mn II 





0.08 



Col 


-0.19 

-0.04 


0.28 

0.22 

0.37 

Cul 



-0.61 

-0.28 




Zn I 

0.23 

0.12 

0.21 

0.29 

0.09 

0.25 

0.22 

Sr II 

-0.45 

-0.78 

-0.35 

-0.19 

-0.85 

-0.57 

-0.41 

YII 

-0.30 

-0.45 

-0.36 

-0.32 

-0.48 

-0.41 

-0.30 

Zr II 

-0.11 

-0.13 

-0.04 

0.04 

-0.26 

-0.15 

0.01 

Ba II 

-0.26 

-0.31 

-0.18 

-0.10 

-0.44 

-0.38 

-0.22 

La II 

-0.06 

-0.18 

-0.16 

-0.15 

-0.29 

0.13 

-0.05 

Nd II 

0.20 

0.07 

0.07 

0.31 

0.05 

0.11 

0.23 

Eu II 

0.26 

0.13 


0.32 

0.07 

0.15 

0.47 

Dyll 

0.45 

0.40 

0.43 

0.69 

0.31 

0.40 

0.33 

Er II 

0.22 

0.08 


0.46 

0.31 

0.22 

0.43 

Yb II 

0.22 

0.12 


0.11 

0.07 


0.15 


^ These abundances are [X/H], instead of [X/Fe]. 


Note. — The machine-readable version of the entire table is available in the online 
journal. A portion of the table is shown here for guidance concerning form and content. 
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Table 9. Average abundances for the individual evolutionary stages and the complete cluster 


M68 


Ion 

RGB 

([X/Fe]> 

RHB 

([X/Fe]> 

BHB 

([X/Fe]> 

M68 

([X/Fe]) Cr([x/Fe]) 

Fe It 

-2.51 

-2.42 

-2.25 

-2.423 

0.141 

Fe lit 

-2.45 

-2.42 

-2.25 

-2.402 

0.156 

OI 

0.69 

0.52 

0.31 

0.506 

0.275 

Nal 

0.20 

0.06 

0.35 

0.165 

0.264 

Mg I 

0.36 

0.35 

0.35 

0.353 

0.170 

All 

-0.23 

-0.40 

-0.07 

-0.276 

0.369 

Si I 

0.40 

0.16 

-0.13 

0.207 

0.238 

Cal 

0.36 

0.43 

0.24 

0.360 

0.111 

Ti I 

0.08 

0.43 


0.215 

0.191 

Ti II 

0.42 

0.43 

0.49 

0.438 

0.134 

Sc II 

0.21 

0.37 

0.22 

0.270 

0.134 

Cr I 

-0.08 

0.19 

0.18 

0.028 

0.165 

Cr II 

0.36 

0.44 

0.33 

0.384 

0.111 

Mnl 

-0.35 

-0.40 


-0.366 

0.150 

Mn II 

-0.11 

0.05 


0.015 

0.340 

Fe I 

0.01 

0.00 

0.10 

0.027 

0.105 

Fell 

0.06 

0.00 

0.10 

0.049 

0.137 

Col 

0.07 

0.25 


0.146 

0.188 

Cul 

-0.54 



-0.540 

0.181 

Znl 

0.19 

0.45 


0.226 

0.133 

Sr II 

-0.56 

-0.16 

-0.78 

-0.444 

0.296 

YII 

-0.42 

-0.29 


-0.363 

0.110 

Zr II 

-0.10 

0.19 


0.032 

0.198 

Ba II 

-0.28 

-0.17 

-0.55 

-0.287 

0.163 

La II 

-0.13 

0.23 


-0.098 

0.164 

Nd II 

0.13 

0.37 


0.166 

0.132 

Eu II 

0.19 

0.29 


0.227 

0.141 

Dy II 

0.38 

0.65 


0.485 

0.215 

Er II 

0.35 

0.77 


0.523 

0.297 

Yb II 

0.15 

0.41 


0.296 

0.173 


^These abundances are [X/H], instead of [X/Fe]. 
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Table 10. Abundance uncertainties 



472 (RGB tip ) 

172 ( RGB ) 

36 ( RHB ) 

337 ( BHB ) 

Species 

a 

a 

G 

G 

EW species 

O I 

± 0.08 

± 0.08 

± 0.14 

± 0.19 

Na I 

± 0.18 

± 0.13 

± 0.20 

± 0.11 

Mg I 

± 0.27 

± 0.17 

± 0.09 

± 0.29 

A1 I 

± 0.11 

± 0.18 

± 0.09 

± 0.13 

Si I 

± 0.15 

± 0.21 

± 0.15 

± 0.09 

Ca I 

± 0.12 

± 0.10 

± 0.11 

± 0.12 

Sc II 

± 0.09 

± 0.18 

± 0.13 

± 0.07 

Ti I 

± 0.16 

± 0.15 

± 0.10 


Ti II 

± 0.08 

± 0.08 

± 0.11 

± 0.11 

Cr I 

± 0.10 

± 0.08 

± 0.20 

± 0.15 

Cr II 

± 0.05 

± 0.08 

± 0.11 

± 0.09 

Mn I 

± 0.10 

± 0.23 

± 0.08 


Mn II 


± 0.24 

± 0.25 


Fe I 

± 0.12 

± 0.13 

± 0.13 

± 0.18 

Fe II 

± 0.16 

± 0.14 

± 0.10 

± 0.13 

Co I 

± 0.18 

± 0.26 

± 0.13 


Cu I 

± 0.08 




Zn I 

± 0.07 

± 0.09 

± 0.13 




Synth species 



Sr II 

± 0.15 

± 0.14 

± 0.13 

± 0.13 

Y II 

± 0.11 

± 0.16 

± 0.17 


Zr II 

± 0.09 

± 0.12 

± 0.09 


Ba II 

± 0.10 

± 0.08 

± 0.13 

± 0.13 

La II 

± 0.10 

± 0.09 

± 0.18 


Nd II 

± 0.16 

± 0.12 

± 0.17 

± 0.16 

Eu II 

± 0.09 

± 0.09 

± 0.20 


Dy II 

± 0.14 

± 0.12 

± 0.17 


Er II 

± 0.16 

± 0.09 

± 0.14 


Yb II 

± 0.14 

± 0.10 

± 0.10 

± 0.17 








